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We present a theoretical study of a possibility of superconductivity in a three dimensional molec- 
ular conductor in which the interaction between electrons in doubly degenerate molecular orbitals 
and an intramolecular vibration mode is large enough to lead to the formation of E ® /3 Jahn- 
Teller small polarons. We argue that the effective polaron-polaron interaction can be attractive for 
material parameters realizable in molecular conductors. This interaction is the source of supercon- 
ductivity in our model. On analyzing superconducting instability in the weak and strong coupling 
regimes of this attractive interaction, we find that superconducting transition temperatures up to 
100 K are achievable in molecular conductors within this mechanism. We also find, for two particles 
per molecular site, a novel Mott insulating state in which a polaron singlet occupies one of the 
doubly degenerate orbitals on each site. Relevance of this study in the search for new molecular 
superconductors is pointed out. 
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I. INTRODUCTION 

Ever since the experimental discovery 1 of supercon- 
ductivity by Jerome et at, in 1980 in a Bechgaard 2 salt 
(TMTTF) 2 PF 6 (with a T c of 1.2 K under a pressure 
of 6.5 Kbar), search for new higher T c molecular su- 
perconductors has been a vigorous field of research 3 . A 
large number of molecular superconductors 3 with vary- 
ing degrees of dimensionality have been discovered since 
1980, and a continuing multidisciplinary search for new 
higher T c molecular superconductors is currently under- 
way. Theoretical studies of a possibility of superconduc- 
tivity in molecular conductors and conducting polymers 
have had a positive influence in the development of this 
field. Indeed, Little's 4 prediction of high superconduct- 
ing transition temperatures within an exciton mediated 
mechanism of superconductivity in a hypothetical con- 
ducting polymer chain (with polarizable molecules peri- 
odically attached to the spine) played a boosting role for 
this field 5 . For a discussion on the future prospects for 
superconductivity in conducting polymers see Heeger 6 . 

One of the strategies in the search for higher T c 
molecular superconductors continues to be to search for 
molecules with large electron- intramolecular vibration 
coupling (EMVC) and to crystallize them with the hope 
that the solid will be metallic or can be made metal- 
lic by charge transfer from another molecule (in charge- 
transfer salts) and/or by application of pressure. If the 
material is metallic and if EMVC is the pairing glue, 
large EMVC implies a large T c within the frame-work 
of BCS theory 7 . In a molecular conductor, the phonon 
spectrum forms two distinct groups: high-frequency in- 
tramolecular phonons and intemiolecular phonons which 
have relatively low frequencies (< 300 cm -1 in the cur- 
rently known organic superconductors). Where calcu- 
lations are available 8 , the EMVC's are larger than the 
electron- intemiolecular vibration coupling, though there 
are exceptions to this trend 9 . But, when EMVC's be- 



come large, there is a possibility of polaron formation 10 , 
a possibility not included in the Migdal-Eliashberg 11,12 
extension of the BCS theory. If EMVC is large enough to 
lead to small polaron formation, a possibility for super- 
conductivity to arise is by the polaron pairing through 
the exchange of low frequency intemiolecular phonons. 
Recently we investigated 14 this possibility in a simple 
model molecular conductor. In that work we considered 
molecular orbitals, overlaps of which on nearest neighbor 
(NN) molecules produces the conduction band, to have 
no orbital degeneracy. 



Now, depending on the symmetry of the molecule, the 
molecular orbitals which participate in the band forma- 
tion can be degenerate or non-degenerate. The purpose 
of the present paper is to investigate the possibility of su- 
perconductivity in a three dimensional molecular conduc- 
tor in which the molecular orbitals participating in the 
conduction band formation are doubly degenerate and 
EMVC between electrons and a single intramolecular vi- 
bration mode is large enough to lead to the formation 
of Jahn- Teller (JT) small polarons. This is the E ® (3 
JT model and is the simplest JT system. Molecules with 
tetragonal symmetry arc examples of such a system. 



The rest of the paper is organized as follows. In Sec. 
II, we introduce our model for the molecular conductor, 
perform a Lang-Firsov transformation on it to obtain the 
JT polaron Hamiltonian. In Sec. Ill, a BCS-mean-field 
(BCS-MF) theory is presented which is applicable in the 
weak coupling regime of the effective attractive interac- 
tion between the JT small polarons. In Sec. IV, we 
analyze the strong coupling Bose regime. An insulating 
state obtained for the case of two electrons per site is 
discussed in Sec. V. The conclusions are given in Sec. 
VI. 
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II. JAHN-TELLER SMALL POLARONS AND 
THEIR EFFECTIVE INTERACTIONS 

We employ the following model to study the E® (3 JT 
polaron formation and their effective interactions: 



* = *EE4 
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In the above 7?, i and j are the molecular site indices, a 
(=1, 2) is the orbital index of the two degenerate molec- 
ular orbitals, t (< 0) the hopping energy between similar 
orbitals on NN molecules, u the molecular vibration fre- 
quency, G the electron-phonon (el-ph) interaction energy, 
M the mass, U the Coulomb interaction strength, and N 
the number of sites. Furthermore, Cj aa and 6 q are the 
destruction operators of electrons and phonons respec- 
tively, and 7ij a(7 — Cj aa Cj a a- A chemical potential will 
be introduced later when we have to fix the electron den- 
sity. In the Appendix, using an argument given in Ref. 15 
we show that, for 2G/w ^> 1, the E®e JT model reduces 
to E ® (3 JT model. A detailed numerical study of po- 
laron and bipolaron formation in E ® e JT model in one 
dimension and for the case of one and two electrons was 
recently published by Shawish et.al. 16 . 

Our calculation proceeds along the following lines. 
First we apply a multi-band Lang-Firsov (LF) 
transformation 17 which produces the JT small polarons, 
reduces the infraorbital Coulomb repulsion, and increases 
the mterorbital Coulomb repulsion. We will argue that 
the intraorbital polaron-polaron interaction can be at- 
tractive for realistic range of frequencies and el-ph cou- 
pling realizable in molecular conductors. This inter- 
action then is the source of superconductivity in our 
model. In the next step, we project out mierorbital dou- 
ble occupancies employing a Gutzwiller approximation 18 
method in the weak and strong coupling regimes of the 
attractive interaction to obtain effective Hamiltonians 
in these regimes. These effective JT Hamiltonians are 
used to study superconductivity in the weak-coupling 
BCS-MF and strong-coupling (Bose) regimes. The su- 
perconducting transition temperatures obtained are then 
analyzed as a function of phonon frequency, el-ph inter- 
action strength, Coulomb repulsion strength, original 
band- width, and doping. We will also see that, for the 
case of two polarons per site, we obtain a novel insulating 
state with a polaron singlet occupying the lattice sites. 
For this paired polaron Mott insulator, we show that one 
obtains a novel Orbital Density Wave State (ODW) for 
large and finite mierorbital repulsion. We also point out 
a mapping of this ODW sate to the Ising model. From 
our analysis of the superconducting transition temper- 
ature in the BCS-MF regime and the Bose regime, it 



will be shown that superconducting T c up to 100 K is 
achievable for this mechanism. Now we go to the details 
of the calculations. 

As stated, we apply a multi-band LF transformation 
to the H. The transformation e s He~ s = Ht transforms 
H into the polaron representation. LF is a unitary trans- 
formation with S = S\ + S2, where S a (a = 1, 2) is 

S a = (-l)^^^n jaff e^'(6 q -&L q ). (2) 



(1) Using this S, we have 



{n ]la - n j2a )e 1 ^ , (3) 
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and 

Xj a = exp 



(5) 



Using the above operators, the transformed Hamiltonian 
is 



where 



Ht — Hke + Hr + Ha + Hj 
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In the above equations 

r 2 

U R = U + 2— , 
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The LF transformation produces three effects. The first 
is the polaronic band-width reduction, of course. The 
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second and third effects are the enhancement of the in- 
ferorbital Coulomb repulsion and the reduction of the 
infraorbital Coulomb repulsion, respectively. These sec- 
ond and third effects have different origins: while the 
second has its origin in the JT nature of the system, the 
third is a small polaron effect. On phonon averaging 19 of 
Ht as usual, neglecting a constant coming from Hi, one 
obtains 

h p ke + h r + h a , (13) 



where 



H 



KE 



in which 



tp = t x exp 



— ) coth ( — 
u I 2 



(14) 



(15) 



and (3 = l/k B T. 

Now, Hp is the Hamiltonian for a collection of JT small 
polarons and these polarons interact among themselves 
though in the intra and inter orbital interactions Ua and 
Up, respectively. At this stage, we demand that Ua be 
an attractive interaction. This is very much realizable for 
the values of U, G, and w in a range possible in molec- 
ular conductors. Some idea about the range of U, G, 
and io would be useful to get a better feel. The high- 
est infromolecular vibration frequency is unlikely to ex- 
ceed the frequency (4161 cm^ 1 ) of the lightest molecule 
(H 2 , the Hydrogen molecule). Among the calculated 
values 20 of G, the highest seems to be around 200 meV 
(obtained for a mode of frequency 1656 cm -1 in Ben- 
zene). Practically all the existing molecular supercon- 
ductors are built from the molecules TMTSF, BEDT- 
TTF, DMIT, DMET, BETS, or C 60 . This is arguably a 
small subset of already synthesizablc molecules. Then, it 
is quite possible that higher values of G are realized in 
many molecules. As for U , large and highly polarizable 
molecules can be expected to have small C/'s. In the exist- 
ing molecular superconductors 3 , U seems to be in a range 
around 1 eV. We note that the U is not for an isolated 
molecule, but for a molecule sitting on a lattice site in 
the solid so that U includes the effects of polarization of 
molecules surrounding a given molecule and consequently 
is reduced from the U value for the isolated molecule. 
Now, if we tentatively fix an upper limit for G around 
600 meV and keep the upper limit of u) as 4161 cm -1 , 
it is not unreasonable to expect C/a to be negative in a 
range of U which obtains in molecular conductors. Given 
the estimates of el-ph coupling to individual modes, this 
situation is very unlikely to be realized in the existing 
organic superconductors. But, that does not preclude its 
realizability in another important class of materials: the 
molecular conductors yet to be discovered. In our study, 
we assume that Ua is negative. This attractive Ua is 
then the origin of superconductivity in our theory. In 
the next two sections, we study superconductivity in the 
weak and strong coupling regimes of Ua- 



III. BCS MEAN-FIELD THEORY IN THE WEAK 
COUPLING REGIME OF U A 



The weak coupling regime of Ua is the range of Ua 
in which the superconducting T c and the pair formation 
temperature are the same. In this section, we project out 
inferorbital double occupancies, using the Gutzwiller ap- 
proximation method, to obtain an effective Hamiltonian 
in the weak coupling regime of Ua and then study the su- 
perconductivity using a BCS-MF theory. The Gutzwiller 
projection renormalizes tp and Ua, and decouples the 
orbitals. First, we recall that the bare el-el repulsions 
[in Eq. (1)] are strongly renormalized [Eqs. (11) and 

(12) ]: while the inferorbital repulsions have increased, 
the infraorbital repulsion is suppressed to an extent that 
it has become attractive (see the previous paragraph). 
The effect of Up is to strongly reduce the probability of 
polarons occupying different orbitals on a given site. On 
the other hand, Ua (now attractive) increases the proba- 
bility of two polarons (with opposite spins) occupying the 
same orbital on a given site. Both these aspects have to 
be considered while doing the Gutzwiller projection. We 
take, for simplicity, the Up to be infinity and completely 
project out the inferorbital double occupancies. Qualita- 
tively one can see that projecting the inferorbital double 
occupancies reduces the polaron band-width (through a 
density dependent function) since it is energetically un- 
favorable for a polaron to hop from an orbital (say 1) 
on a site to a target site where the orbital 2 is occupied. 
Also, strong reduction of inferorbital double occupancies 
increases the infraorbital double occupancies and conse- 
quently modifies the infraorbital attraction. We do the 
inferorbital double occupancy projection in an approxi- 
mate way. First consider the part H^ E + Hp of Hp [Eq. 

(13) ]. We write the Gutzwiller wave function \ipGw) as: 



\ipGw) = - (1 - v)nnan i2 a'}\ipo) , 



(16) 



where 



nm= n 4ii«>® n c i 2 io) 



(17) 



k 2 <k F 



For calculating the expectation value of H^ E and Hp 
in \tpGw), we use a simple intuitive method proposed 
by Okabe 21 . Since the maximum number of polarons 
(Sao- n au) a site can accommodate is four, there are six- 
teen possibilities for occupancy for finite Up. As men- 
tioned earlier, we take Up infinity limit for simplicity 
[t) = in Eq. (16)], and consequently the sixteen pos- 
sibilities reduces to seven. In the paramagnetic case 
(i.e., ni| = nil = n 2| = n 2i = n/4), the probability 
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for a site to be vacant is [1 — (n — 2d)] and the probabil- 
ity for a site to be singly occupied is (n/4 — d). Here d 
is the infraorbital double occupancy. Then, considering 
the hopping process given in Fig. 1, the Gutzwiller band 
narrowing factor is 



q(n,d) = WWEME^KlVMES , (18) 

4 V 1 4 > 

where the numerator and the denominator corresponds to 
hopping in the presence and absence of mf erorbital repul- 
sion, respectively. The ground state energy for H^ E +Hn 



is 



E c (n,d) = q(n,d)^2< 



(19) 



where e arT [= e] is the average kinetic energy per polaron, 
and is independent of the indices a and a. Minimiza- 
tion of E c (n,d) with respect to d leads to the equation 
to determine d for a given n: 



(4 - 6n)d 2 + (1 - 2n + \n 2 )d + 



n 3 - n 2 
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0. 

(20) 



For n — 2, the above equation implies that d = 0.5, which 
means that one of the orbitals on all the sites are occu- 
pied by a singlet. For a general value of n, one has to 
solve Eq. (20) to obtain the variation of d as a function 
of n. The optimized values of the band narrowing factor 
q(n, d op t) and the infraorbital double occupancy {dopt) as 
a function of n are shown in Fig. 2 and Fig. 3, respec- 
tively. From the numerical results shown in the figures, 
we have 



q(n) = q(n,d opt ) 



and 



d, 



opt 



(21) 



(22) 



Furthermore, it is clear that as n — > 2, the polaron band- 
widths are narrowed and a Mott insulator is obtained for 
n = 2. This Mott insulator has a singlet in one of the or- 
bitals on each lattice site. The properties and excitations 
of this insulator will be discussed separately in Sect. V. 
Next consider the term H A in Eq. (13). To get the total 
ground state energy, we add to E c (n, d) the contribution 
from H A obtained using \ipGw) calculated for optimized 
d (d op t). This contribution per site is 

-\Ua\ 7T T7 r =-2\U A \dopt- (23) 

VPGWWGW) 

So, the total minimized energy per site is 



E (n) = ^2 l( n ' d o P t)e a a ~ 2\U A \d op t ■ 



(24) 



Using Eqs. (21) and (22), the minimized energy per site 
is 



E(n) w2(2-n)e- 



Wa\i 



(25) 



Using the results of the Gutzwiller projection of inter- 
orbital double occupancies, we write an effective Hamil- 
tonian for the JT polarons as: 



Hp = q(n)t P J2Y, c i 



UA^2nj a ^nj a i . (26) 



The effects of Gutzwiller projection are the density de- 
pendent band-width renormalization and a renormaliza- 
tion of the on-site infraorbital attractive interaction from 
U A to U A . From Eq. (25), it follows that U A = 2U A . We 
notice that this enhancement is similar to the enhance- 
ment of antiferromagnetic exchange interaction (J) ob- 
tained on Gutzwiller projection in a t— U — J model stud- 
ied in Ref 22 . For the t — U — J model, they found that the 
bare J steadily increases with increasing electron density 
to reach 4J in the Mott insulator at the half-filling of 
their non-degenerate band. Notice also, for our model, 
that orbitals are decoupled on Gutzwiller projection. 

Now, H is clearly an effective attractive Hubbard 
model. We note that, unlike the usual phonon medi- 
ated attractive interaction which is attractive only in the 
Debye shell around the Fermi surface, all the polarons 
experience the attraction U A or U A . In the next section, 
we study superconductivity in the weak coupling regime 
of U A using a BCS type mean-field (BCS-MF) analysis 
to obtain superconducting energy gap and T c . 

To move forward, we have to fix the weak and strong 
coupling regimes of the single band attractive Hubbard 
model (H AHM = tj^ijo- c L c j<? ~ l^l Z)j n rt n jl) in terms 
of |E/|/t. We can approximately fix these ranges us- 
ing the existing studies 23-25 on three dimensional non- 
degenerate attractive Hubbard model. A comparison 
of results for superconducting T c obtained using three 
methods (T-Matrix Approximation, Determinant Quan- 
tum Monte Carlo, and Dynamical Mean Field Theory) 
and the BCS-MF is given in Fig. 2 of Ref. 25 Notice that 
these results are for quarter filling of the band. The 
regime of weak coupling, for a given n, is the range 
of \U\/t within which the pair formation temperature 
and the superconducting transition temperature are the 
same. Using the calculation of Pauli susceptibility (which 
would be suppressed above the superconducting T c in the 
case pairs form above T c ), the authors of 25 argued that 
< \U\/t < 4 is the BCS-MF regime for n = 0.5. Con- 
sidering together all these results, we tentatively fixed 
the BCS-MF regime of the 3D attractive Hubbard model 
to be < \U\/t < 3. As for the strong coupling regime, 
all the electrons are in paired states below a pair forma- 
tion temperature of 0(|C/|). We fixed the Bose regime 
to be for \U\/t > 12. The intermediate coupling regime 
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then is obviously 3 < \U\/t < 12. In the weak coupling 
regime of Ua, the ratio equivalent to \U\/t in our case 
is, from Eq. (26), \U a\ I q{n)t p . In the strong coupling 
regime of Ua studied in the next section, the Gutwiller 
projection has to be done in a different way than that 
given in this section, and that gives \Ua\/ \Zq{n)tp [where 
q(n) = (1 — n/2)/(l — n/4)] as the ratio equivalent to 
\U\/t. We can now proceed to the analysis of supercon- 
ductivity in the BCS-MF regime of the effective attrac- 
tive interaction Ua- 

We have from Eq. (26), 



(27) 



where 



and 



&(k) = q(n)e P (k)-», 



(28) 



e P (k) = 2tp[cos(k x a) + cos(k y a) + cos(k z a)] , (29) 

for a simple cubic lattice. We have introduced a chem- 
ical potential (fi) to fix the number density. Given the 
form of Ua, only isotropic s-wave pairing is possible. A 
BCS-MF theory of Hf, gives: 

H%=Y, &( k )*W + ^(A* C_ kQi C kQT + h.C.) , (30) 
kcr k 

where the superconducting gap parameter is 



k = uaJ2(cLA 



-kal 



(31) 



The superconducting gap and the particle density is de- 
termined from 



U A 



+q(n)D P jVp( e a ) 
-q(n)Dp Q(n) 



F(e a ,A a ,^(3)de a , (32) 



and 
1 



f +q(n)D P N , a \ 

J- q {n)D P Q{n) 



(33) 



where 



tanh (.5/V(e a -M) 2 + |A Q | 2 ) 
F(e a ,A Q , M ,/3) = V - >-. (34) 

Here A P (e Q ) is the polaron Density Of States (DOS) 
corresponding to the polaron band ep(fc). Assuming 
a square one-polaron single-spin DOS (i.e., Np(e a ) — 
l/2Dp = l/12tp), the solutions of the gap and chemical 
potential equations give: 



A a (T = 0) 



1 T6 p 

sinh ( j) 



and 



k B T c 
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1.13-Dpe - * 



(35) 



(36) 



In the above equations, A = 2\U — 2g 2 uj\/ (2q(n)Dp), 
Dp = D Q e~ 9 , 5 = 2 — n, g = G/co, and n is the number 
of polarons per site. For the values of g and uj we use 
in numerical calculations, the temperature dependence of 
the polaron band-width is practically nil for the the tem- 
perature range (< 300 K) of interest to us. The equation 
for T c is valid for [(q(n)D P ± fi)/k B T c ] > 1. See also the 
note in Rcf 26 . 

Now we go to the estimates of A a (0) and Tc con- 
sidering realistic values of g, uj, D 0l and U applicable 
to molecular conductors. While making the estimates, 
one has to satisfy several constraints. These are: (i) JT 
small polaron formation condition (g 2 ui > q(n)D ); (ii) 
Phonon frequency should be below the frequency (4161 
cm -1 ) of the lightest molecule (H2); (iii) A realistic up- 
per limit for el-ph interaction strength was earlier fixed at 
600 meV; and (iv) Range of validity of BCS-MF theory 
requires that A < 0.25. These conditions enforce severe 
constraints on the parameters g, uj, D , and U for which 
the system can undergo superconducting instability and 
if that happens on the achievable T c 's. The variation of 
A Q (at T = 0) with 6 is shown in Figs. 4-6, for dif- 
ferent values of uj, D q , and g, respectively. Since the 
trends are similar for A a and T c , we will discuss the 
T c curves shown in Figs. 7-9. The lower limit of <5 for 
the curves in Figs. 4-9 is determined by the condition 
on A, and the upper cut-off as a function of S enforces 
the small polaron condition. The variation of T c with 
S for various values of uj is shown in Fig. 7. For fixed 
U, D and g, BCS-MF superconductivity is confined to 
the uj range: u> min < uj < uj max where uj min = U/2g 2 
and uj max = uj min + (0.25q(n)Dp/2g 2 ). In this allowed 
range of u>, moderate values of T c are possible. Increas- 
ing uj is found to increase the T c 's achievable and expands 
the density range in which superconductivity is possible. 
The variation of T c with <5 for different values of D Q is 
shown in Fig. 8, and it shows that in the range shown 
changes in D Q do not have much effect. On the other 
hand, the T c vs. S curves shown in Fig. 9 show that 
small changes in g severely narrow the range of S for 
which superconductivity is possible. It is clear that the 
various conditions mentioned earlier severely limit the 
range of 6 for which superconductivity is possible and 
the values of T c 's achievable. We would also like to note 
that below the lower cut-off of 5 for each curve in the 
Figs. 4-9, the A moves out of the BCS-MF regime into 
the intermediate to strong coupling regime and whether 
superconductivity is possible in that range of S cannot be 
addressed using the weak coupling formulas used in this 
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section. The values of U, g, D Q , and u> used in making 
the T c estimates are definitely not unrealistic considering 
the variety of molecules synthcsizable and crystallizable. 
Though band-filling control has not yet been achieved in 
molecular conductors, such a possibility cannot be ruled 
out in the future. Even if band-filling control is not pos- 
sible, our results show that JT polarons can undergo su- 
perconducting instability with moderate T c 's in a realistic 
range of U , g, D Q , and 10 realizable in molecular conduc- 
tors. This completes the analysis of superconductivity in 
the BCS-MF regime of our model. In the next section, 
we consider the strong coupling regime (the Bose regime) 
of U A . 



The hopping processes corresponding to the first and sec- 
ond terms in E 2 are shown in Fig. 10. E 2 can be cal- 
culated considering the probability amplitudes for the 
processes shown in the figure and we obtained it to be 
the ground state energy of the Hamiltonian 



ij a 

f 2 pg(n) vv Zt 2 p q(n) ^ 

ij aa jctcr 



IV. BOSE CONDENSATION IN THE STRONG 
COUPLING REGIME OF U A 

In this regime, all the JT polarons are in paired states 
below a temperature of 0(|L/^|). To obtain the Bose 
condensation temperature of these pre- formed pairs, we 
start from the Hamiltonian Hp [see Eq. (13)]. Since 
the largest terms in Hp are Hp and Ha (in the strong 
coupling regime), we treat (Hp + Ha) as the part which 
determines the ground state energy and take H^ E as a 
perturbation. Now, in the ground state of (Hp + H A ), 
a site has either a singlet occupying one of the orbitals 
or is empty. The mferorbital pairs are projected out be- 
cause of the strong mterorbital repulsion (Up). When 
we switch on the hopping term (H^ E ), the pairs become 
mobile through a second order (in H^ E ) process. For 
the ground state of (Hp + H A ) we have 



(H R + H A )\p)=E \p), 



(37) 



where E = (— \U A \/2) X)j Q(T n jaa is the ground state en- 
ergy. The second order pair hopping term then is, 



Ei 



t 2 

Lp 



I [/J ^Vl^^~iaa^ 

1 1 r ij aa 

x EE c iW' c ™^'ip> • 

Ira 0<j' 



r)(r\ 



(38) 



In Eq. (38), the r sum is over states other than the degen- 
erate ground states. Since ^ r \r)(r\ = I — J2 P \p)(p\ an d 
noting that (p\T\p') = 0, the second order term becomes 

t 2 

E 2 = -y^jd'lEE^^EE^V' " 1 / 3 -'^ • 



Im (3 a' 



(39) 

The E 2 is non-zero for (A) I = j; i = m; a = (3; a = a' 
or (B) i = I; j = m; a = (3; u = —a 1 . Then, 



E 2 = 



\Ua\ 

(piE 



((p| 2 E c LT c Li c j«l c j«Tb> 



i] a 

C iaa C jaa c j aa Ci a(J 



If))- 



(40) 



where q(n) = (1 — n/2)/(l — n/4) is the Gutzwillcr pro- 
jection factor and Z is the co-ordination number of a 
site in the lattice. The Bose condensation temperature 
(T B ) in the strong coupling regime of a three dimensional 
non-degenerate attractive Hubbard model was obtained 
in Ref. 27 . Comparing our results with theirs, one can 
immediately write down the Tp for our case, and it is 



k B T B = 



2(n-2)q(n)Zt 2 



(42) 



The above formula is relevant for \U a\ I (\J q(n) tp) > 12. 
Here it must be pointed out that the finite tempera- 
ture analysis in this section was done in real space and 
hence one need not worry about the orthogonality of 
the Gutzwiller projected excited state wave functions ob- 
tained by varying the occupation numbers rikaa, as was 
done in Ref. 26 . Note that for n = 2, we have a Mott 
insulator and Tp is zero. While making Tp estimates, 
we have to again satisfy the constraints (i)-(iii) given in 
the previous section. The small polaron formati on co n- 
dition in the strong coupling regime is g 2 u > ^q(n)D 
where D is the original half-band-width equal to 6t for 
a simple cubic lattice we have considered. The variation 
of T B as a function of S for various values of uo is shown 
in Fig. 11, and it is seen that decreasing ui increases the 
Tb's and contracts the range of 5 in which Bose conden- 
sation is possible. The upper cut-off in Figs. 11-13 is 
due to the violation of the small polaron condition. In 
Fig. 12, we have displayed Tp vs. 5 for different values 
of D Q . Decreasing D decreases Tp while expanding the 
S range for Bose condensation. Finally, Tp vs. 8 for var- 
ious values of g is shown in Fig. 13. Increase of g is seen 
to sharply reduce Tp . Overall, it is clear from the figures 
that Tg's up to 100 K can be possible in the strong cou- 
pling regime. Furthermore, we note that slightly higher 
values of Tp can be obtained in the intermediate cou- 
pling regime as compared to the strong coupling regime 
threshold (|I7|/i w 12) (see Fig. 2 in Ref. 25 ). 
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V. THE MOTT INSULATOR WITH SINGLETS 

We have noted earlier that the ground state of the 
Hamiltonian Hp is a Mott insulator for the case of two 
electrons per site. This Mott insulating state obtains 
independent of the value of Ua (< 0) (see also the 
note in Ref. 28 ). The insulating state is obtained be- 
cause of the strong mierorbital repulsion Ur between 
the JT polarons. The repulsion Ur narrows the po- 
laron band as the density increases and eventually drives 
it zero at n = 2. Now, due to the presence of finite 
infraorbital attraction between the polarons, the insu- 
lating state has a polaron singlet occupying one of the 
orbitals on site. Since the singlet can occupy either of 
the two orbitals available on a site, the insulting ground 
state is highly degenerate when Ur is infinity. When 
Ur is large and finite, a pair on a site can gain energy 
by virtual hopping to a nearest neighbor vacant orbital. 
This gain in energy, in second order perturbation theory, 
is (—Ztp)/(2UR + \Ua\) per polaron. Hence for finite 
and large Ur the degeneracy of the insulating state re- 
duces to two and the singlets order in alternate orbitals 
on NN sites and thus we have an Orbital Density Wave 
state. We also note that one can map the paired po- 
laron Mott insulator to an Ising model since at each site 
the pair could be in cither of the orbitals. Assigning 
pseudo-spins j to occupation of orbital 1 and j to that 
of orbital 2, the Ising model is H = J2<ij> JSi-Sj with 
J = tp/(2UR + \Ua\)- The ground state corresponds to 
the orbitally antiferromagnetic state mentioned earlier. 
When the electron density deviates from n = 2, the in- 
sulator undergoes an insulator to metal transition and 
the superconducting instability of this metallic state was 
analyzed in the previous sections. 



VI. CONCLUSIONS 

In this paper we presented a theoretical study of su- 
perconductivity in a three dimensional molecular conduc- 
tor in which the molecular orbitals participating in the 
conduction band formation are doubly degenerate and 
EMVC for electrons in these orbitals interacting with 
non-degenerate intramolecular vibration mode is large 
enough to lead to the formation of Jahn- Teller ( JT) small 
polarons. We argued that the effective polaron-polaron 
interaction can be attractive for material parameters re- 
alizable in molecular conductors, and this interaction is 
the source of superconductivity in our model. On ana- 
lyzing superconducting instability in the weak and strong 
coupling regimes of this attractive interaction, we found 
that superconducting transition temperatures up to 100 
K are achievable in molecular conductors within this 
mechanism. We also find, for two particles per molec- 
ular site, that the ground state is a novel Mott insulat- 
ing state in which a polaron singlet occupies one of the 
doubly degenerate orbitals on each site. In the infinite 



interorbital repulsion case, this Mott insulator is highly 
degenerate since the singlet has the freedom to occupy 
either of the orbitals on a site. On the other band, in 
the case of large but finite inierorbital repulsion case, we 
found that the degeneracy of the ground state is reduced 
to two with the singlets orbitally antiferromagnetically 
ordered. When the number of particles per site deviates 
from two, the Mott insulator undergoes an insulator to 
metal transition. 
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APPENDIX A 

Here we show that, for strong coupling (i.e., 2G/oo ^S> 
1), E ® e JT model is equivalent to the E <g> (3 model in 
Eq. (1). For E <g> e, 



2 Z-^-*"* 2 

ia k\ 

- G V2M^E(4.4.)^ 0li Q 

ia 

+ U E n i«T n iai + u E n i^ n J^ 



2i \ ( dilcr 
Q2i —Qli J \ di2a 



(Al) 



Now we do a rotation transformation on H . For this 
purpose, let 

Qu Q2i\_ n ( -cos(6i) sin(6i)\ 
Q 2 * -Qu )~™\ sin(6i) cosiOi) ) ' [AZ > 

where Qi = Q 2 U + Q\ and 

r,,„ \ / sin(?f) cos(?f) \ ( dn 



Ci2a J \ COS(^) -Sin(^) J \ di2a 

Using the above, H is transformed to: 



■ (A3) 



*=*e(4„cl)(;; 

ija 



Cjlcr 
Cj2a 



ia 

Muj 2 x - 2 _L_ x - 



t J \ ( Qi \ / Cii 



i2 °J V -Qi M c 42ct 



d 2 I d_ 1 d 2 

2M ^ yoo 2 + q~ oq~ + qI oe 2 

+ Uy^njainjai + E n i^nj2u' ■ 

ja jcrcr' 



(A4) 



As argued in Ref. 15 , for 2G/u > 1, {l/Qi){d/dQi) 
and (1/Qf ){d 2 /d9 2 ) are unimportant, and we are left 
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with an effective single mode equation. Then using 
\j2MujQi = bi + b\, one gets the Eq. (1) given in the 
text. Here it should be pointed out that the single normal 
mode approximation for E®e works only for the low lying 
excitations of the mode. For the transition temperatures 
obtained and the frequencies considered ui/ksTc » 1, 
and hence phonon averaging after the Lang-Firsov trans- 
formation yields for the single mode approximated E <g> e 
system the same result given by Eq. (15) in the text. 



1 D. Jerome, A. Mazaud, M. Ribault, and K. Bechgaard, J. 
Phys. (France) Lett. 41, L95(1980). 

2 K. Bechgaard, C. S. Jacobsen, K. Mortensen, H. J. Peder- 
sen, and N. Thorup, Solid State. Commun. 33, 1119(1980). 

3 T. Ishiguro, K. Yamaji, and G. Saito, Organic Supercon- 
ductors (Springer, Berlin, 1998). 

4 W. A. Little, Phys. Rev. 134, A1416(1964). 

5 D. Jerome, Science 252, 1509(1991). 

6 A. J. Heeger, Rev. Mod. Phys. 73, 681(2001). 

7 J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Phys. Rev. 
108, 1175(1957). 

8 A. J. Berlinsky, J. F. Carolan, and L. Weiler, Solid State 
Commun. 15, 795(1974); N. O. Lipari, M. J. Rice, C. B. 
Duke, R. Bozio, A. Girlando, and C. Pecile, Intl. J. Quan. 
Chem. QM Symp. 11, 583(1977); K. Yamaji, Solid State. 
Commun. 61, 413(1987); J. C. R. Faulhaber, D. K. K. Ko, 
and P. R. Briddon, Syn. Met. 60, 227(1993); V. M. Yart- 
sev and A. Graja, Intl. J. Mod. Phys. B 12, 1643(1998); 
A. Devos and M. Lannoo, Phys. Rev. B 58, 8236(1998); T. 
Kato and T. Yamabe, J. Chem. Phys. 115, 8592(2001); 
T. Kato, K. Yozhizawa, and K. Hirao, J. Chem. Phys. 
116, 3420(2002); J. Sinoa, J. Schliemann, A. S. Nunez, 
and A. H. MacDonald, Phys. Rev. Lett. 87, 226802(2001); 
C. M. Varma, J. Zannen, and K. Raghavachari, Science 
254, 989(1991); V. P. Anthropov, O. Gunnarson, and A. 
I. Liechtenstein, Phys. Rev. B 48, 7651(1993); O. Gunnar- 
son, Rev. Mod. Phys. 69, 575(1997); M. Cote, J. S. Gross- 
mann, M. L. Cohen, and S. G. Louie, Phys. Rev. Lett. 81, 
697(1998). 

9 R. Ramakumar, Y. Tanaka, and K. Yamaji, Phys. Rev. B 
56, 795(1997); R. Ramakumar and K. Yamaji, Eur. Phys. 
J. B 24, 71(2001). 

10 J. Appel, in Solid State Physics, Vol. 21, eds. H. Ehrenre- 
ich, F. Seitz, and D. Turnbull (Academic Press, New York, 
1967), pl93; A. S. Alexandrov and N. F. Mott, Polarons 
and Bipolarons (World Scientific, Singapore, 1995). 

11 A. B. Migdal, Zh. Eksp. Teor. Fiz. 34, 1438(1958) [Sov. 
Phys. JETP 34, 996(1958)]. 

12 G. M. Eliashberg, Zh. Eksp. Teor. Fiz. 38, 966(1960) 
[Sov. Phys. JETP 11, 696(1960)]; 39, 1437(1960)[12, 
1000(1961)]. 

13 T. Holstein, Ann. Phys. (N. Y.) 8, 325(1959); 8, 343(1959). 

14 R. Ramakumar and S. Yarlagadda, Phys. Rev. B67, 
214502(2003). 



15 M. D. Sturge, in Solid State Physics, Vol. 20, eds. H. Ehren- 
reich, F. Seitz, and D. Turnbull (Academic Press, New 
York, 1967), p 91. 

16 S. El Shawish, J. Bonca, Li-Chung Ku, and S. A. Trugman, 
Phys. Rev. B 67, 014301(2003). 

17 I. G. Lang and Yu. A. Firsov, Zh. Eksp. Teor. Fiz. 43, 
1843(1962) [Sov. Phys. JETP 16, 1301(1962)]. 

18 M. C. Gutzwiller, Phys. Rev. Lett. 10, 159(1963); Phys. 
Rev. 134, A923(1964). 

19 S. Yarlagadda, Phys. Rev. B 62, 14828(2000). 

20 T. Kato and T. Yamabe in Ref 8 . 

21 T. Okabe, J. Phys. Soc. of Japan 65, 1056(1996). 

22 F. C. Zhang, C. Gros, T. M. Rice, and H. Shiba, Super- 
cond. Sci. Technol. 1, 36(1988). 

23 R. Micnas, J. Ranninger, and S. Robaszkiewicz, Rev. Mod. 
Phys. 62, 113(1990). 

24 M. Keller, W. Metzner, and U. Schollwock, Phys. Rev. 
Lett. 86, 4612(2001). 

25 A. Sewer, X. Zotos, and H. Beck, Phys. Rev. B 66, 
140504(2002). 

26 The reader may be concerned about the validity of calculat- 
ing T c from a theory on an effective Hamiltonian obtained 
using Gutzwiller approximation. Since BCS-MF T c is much 
smaller than the effective Fermi temperature, the correc- 
tions brought about by the modification required in the oc- 
cupation number (nkaa) to incorporate the finite tempera- 
ture effects in Gutzwiller projection is expected to be small. 
See: T. M. Rice, K. Ueda, H. R. Ott, and H. Rudigier, 
Phys. Rev. B 31, 594(1985); K. Seiler, C. Gros, T. M. 
Rice, K. Ueda, and D. Vollhardt, J. Low. Temp. Phys. 64, 
195(1986). 

27 S. Robaszkiewicz, R. Micnas, and K. A. Chao, Phys. Rev. 
B 23, 1447(1981). 

28 For n = 2, a Mott insulating state is also possible for 
< Ua < Ur. 



_ti_ 



_ti_ 



_ti_ 



_ti_ 



FIG. 1. The inter-site hopping processes considered (see 
text). 
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FIG. 2. The variation of the band-width renormalization 
factor as a function of particle density. 
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FIG. 4. Aq, vs. 5 for various values of w (in eV): 0.348 
(left), 0.352 (middle), 0.355 (right). Values of the other pa- 
rameters are: U = leV, D a = 0.6 eV, and g — 1.2. 




FIG. 5. A a vs. S for various values of D (in eV): 0.70 
FIG. 3. The variation of the infraorbital double occupancy (left), 0.65 (middle), 0.60 (right). Values of the other param- 
as a function of particle density. eters are: (7 = 1 eV, u> = 0.355 eV , and g = 1.2. 
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FIG. 6. A a vs. 8 for various values of g: 1.200 (left), 1.202 
(middle), 1.204 (right). Values of the other parameters are 
(in eV): U = 1, w = 0.355, and D D = 0.6. 
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FIG. 8. T c vs. 8 for various values of D a (in eV): 0.70 
(left), 0.65 (middle), 0.60 (right). Values of the other param- 
eters are: U = 1 eV, u> = 0.355 eV, and g = 1.2. 
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FIG. 7. T c vs. 8 for various values of w (in eV): 0.348 (left), 
0.352 (middle), 0.355 (right). Values of the other parameters 
are: U = 1 eV, D a = 0.6 eV, and g = 1.2. 
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FIG. 9. T c vs. 8 for various values of g: 1.200 (left), 1.202 
(middle), 1.204 (right). Values of the other parameters are 
(in eV): U = 1, w = 0.355, and D = 0.6. 
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FIG. 12. T s vs. 8 for various values of D a (in eV): 0.60 
(top), 0.55 (middle), 0.50 (bottom). Values of the other pa- 
rameters are: U = 0.9 eV, u> = 0.4 eV, and g = 1.2. 
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FIG. 11. Ts vs. 8 for various values of w (in eV): 0.350 
(left), 0.375 (middle), 0.400 (right). Values of the other pa- 
rameters are: U = 0.9 eV, D = 0.6 eV, and g — 1.2. 
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FIG. 13. Ts vs. 8 for various values of g: 1.20 (top), 1.3 
(middle), 1.4 (bottom). Values of the other parameters are (in 
eV): U = 0.9 (thick lines) and U = 1.0 (thin lines), w = 0.4, 
and Do = 0.6. 
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